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Abstract 
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*~~5 The Skyrme model is a non-linear field theory whose solitonic solutions, once quan- 

ta tised, describe atomic nuclei. The classical static soliton solutions, so-called Skyrmions, 

have interesting discrete symmetries and can only be calculated numerically. Mathe- 
i— i matically, these Skyrmions can be viewed as maps between to two three-manifolds and, 

as such, their stable singularities can only be folds, cusps and swallowtails. Physically, 
the occurrence of singularities is related to negative baryon density. In this paper, we 
D calculate the charge three Skyrmion to a high resolution in order to examine its singu- 



larity structure in detail. Thereby, we explore regions of negative baryon density. We 
also discuss how the negative baryon density depends on the pion mass. 



> 

0\ 1 Introduction 

The Skyrme model is a (3 + l)-dimensional nonlinear theory of pions and was conjectured by 
Skyrme as a model of baryons [lj . Subsequently, Witten [2] derived it as an effective action 
of QCD - in the large colour limit. The minimum-energy solutions of the model are called 
Skyrmions and are identified as baryons. In this article we are interested in static solutions 
. ^ of the Skyrme model, best defined by the static energy functional 

& E= 1^J {-^ Tr (^)-T^ Tr ([^^i][^'^']) + m2Tr ( J 2-^)} rf3a; ' ( L1 ) 

where is an sii(2)-valued current, = (8^11)11^ and U(pc,t) is an SU(2)-valued scalar 
field. The parameter m is related to the physical pion mass via m = 2m 7r /(e-F 7r ) where 
F n is the pion decay constant and e is the Skyrme constant. Here we have presented the 
energy functional in so-called Skyrme-units, see [31 H] for the "standard values" and [HI [6] for 
a more detailed discussion of our current understanding. For non-zero pion mass the finite 
energy requirement forces the field U : 1R 3 — > SU(2) to satisfy 

lim U = h, 

r— >oo 
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and we choose the same boundary condition for m n = 0. This one-point compactifies M. 3 to 
the 3-sphere S 3 . Using the fact that the group manifold of SU(2) is also S 3 the field can be 
extended to a map U : S 3 — > S 3 . All finite energy field configurations Z7(x) belong to an 
element of ^(S" 3 ) = Z, and hence have an associated integer B G Z. It is this B which is 
identified as the baryon number and can be explicitly calculated as 



B 



l — J eijkTriRiRjRjdPx, (1.2) 
det(J(x))d 3 a:, 



24vr 2 
1 

2^2 

where J(x) is the Jacobian of the map [7J. Naively, the density of B could be assumed to 
be positive for all x e M. 3 , but this does not have to be the case. It was shown in [8] that 
as 1R 3 covers S 3 - the group manifold of SU(2) - there is a folding structure. On these folds 
the Jacobian-determinant becomes zero which results in tubes of zero baryon density. Inside 
these tubes singularity theory predicts regions of negative baryon density, but this was never 
observed in numerical calculations. This is what we are interested in this paper. It is useful 
to define two quantities, 

B+(x) = ^max{det(J(x)),0}, (1.3) 

8-0*0 = — max{-det(J(x)),0}, 
2tt z 



where B = J B(x)d 3 x and trivially <B(x) = <6+(x) — <B_(x). We refer to these quantities as 
the positive baryon density and negative baryon density, respectively. 

The paper is laid out as follows. Section [2] reviews the geometric formulation of the Skyrme 
model proposed in [7J and discusses the Jacobian. In particular, we show an important 
identity relating the Jacobian of the Skyrme field to a simpler quantity. Section [3] reviews 
the rational map ansatz [9] as well as the non-holomorphic rational map ansatz [8] . In section 
[4] we examine the folding structure of the B = 3 Skyrmion using both full field simulation 
and the non-holomorphic rational map ansatz. Section 4.1 discusses the effects of the pion 
mass term on the folding structure. In section 4.48 we locally expand the pion fields around 
the origin and reproduce the singular surface. We end with a conclusion. 



2 Geometric setting for Skyrmions 

Instead of considering the Skyrme model in a physical field theoretic setting, it can be very 
insightful to consider it geometrically [7J [10]. In the following, we give a brief account of this 
approach and set up our notation. 

A field configuration is a map it from a physical space M. 3 to a target space SU(2). 
Physical space M. 3 and target space SU{2) = S 3 both are 3-dimensional, connected and 
orientable Riemannian manifolds. Here we choose the Kronecker delta 8ij as the flat-metric 
on M. 3 and 7 a/ g as the metric on S 3 . We denote a point in M. 3 as x and a point in S 3 , the 
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image of x, as tt(x). As M. 3 is flat we trivially choose its dreibein to be 8ij, and we choose 
£ M a (7r(x)) as the dreibein on S 3 with 

CfCS-rap = <W (2-4) 

Note the inverse of C/ is C a , i.e. CaCi/* = and CaC/ = S a p . 

Now we can define the Jacobian matrix associated with the map tt(x) as 

J/(x) = (^vr"(x))C^(7T(x)). (2.5) 

The matrix J m M (x) is a measure of the deformation induced by the map 7r at the point x in 
1R 3 . We can define a useful quantity called the strain tensor Z)y, as 

Ai(x) = J>(x)J/(x)5^ = (5 i vr a (x))(^7r' 3 (x))7^. (2.6) 

The strain tensor is invariant under rotations in target space (i.e. rotations of the frame 
fields C M a (7r(x))), but not under orthogonal rotations of the physical space 1R 3 . But it is well 
known that the characteristic polynomial P = det(D — A/3) is invariant under rotations. So, 
we can now define the three invariants, 

Tr(-D) = A 2 + A 2 + A 2 , (2.7) 

\{T,{D)f -\tt{D 2 ) = A 2 A 2 + A 2 A 2 + A 2 A 2 , 

det(D) = A 2 A 2 A 2 , 

where X\ , X\, X\ are the non-negative eigenvalues of the symmetric matrix Dij. 

Reparametrising the SU(2)-valued field in the traditional way with the three Pauli ma- 
trices, r a , and the four scalar fields cx(x),7r a (x) (a = 1,2,3) as 

C/(x) = o-(x) + mt(x) • r, (2.8) 

we see that cx(x) and 7r(x) must satisfy the constraint cr 2 + 7r • 7r = 1 for all x e M. 3 . As 
pointed out in [TU] 7r a (x) in the above geometric discussion can be identified with the vector 
ir a in (2.8) and cr(x) is a function of 7r a (x) which ensures that U(x) G SU(2). Now, we can 
define the induced-metric on S 3 as 

7 Q/3 (x) = <J a/3 H — . 

A short calculation in 1 10 shows that 

Ai(x) = ~Tr(/^), 



and the energy functional (1.1) can be rewritten as 



E = Tr(D) + ^(Tr(D)) 2 — ^Tt(D 2 ), (2.9) 

— ^1 + ^2 + ^1 + ^1^2 + ^1^3 + ^2^3- 



Also, it is easy to see from the above relations that the baryon number integral (1.2) can be 
written in terms of these eigenvalues as 

B = 2^ I AiA 2 A 3 d 3 x. (2.10) 
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2.1 The Jacobian 

In the following we calculate the Jacobian in stereographic coordinates. This is needed for 
the analysis of the numerically found solutions. Consider the Skyrme field in sigma model 
coordinates (cr(x), 7Tj(x)) with a 2 +ir 2 = 1. Then we can define stereographic coordinates by 
projecting from the North pole iV (er = 1) as 

n = T^- (2.H) 

1 — a 

Note that the metric in stereographic coordinates for this chart is given by 

3 . 

^ = E /up 2^ 2 W) 2 , (2-12) 



where = (®n) = iir- Since the metric (2.12) is diagonal we can define the frame 

a=l 

fields 2 

CV = 1 + * N and its inverse ( N » a = 2 ^ 2 5» a . (2.13) 

l 1 + it at 

Hence, the Jacobian is given by 



d$% . „ « r 

2 • 



JNm 11 - JT^-CN^a ~ JT^~ 2 - ( 2 -14) 

OX m OX m 1 + K]\f 



Using (2.6) we can calculate the strain tensor D mn and obtain 

DNmn -(l + R N *)*d^-d^ 5ap - (2 - 15) 

Note W is not an index to be summed over. 

To obtain a well-defined 5*0(3) frame bundle the frame fields (s have to be chosen as 

1 -\- R 2 

(s^ a = ^ ^ a an< ^ ^ s i nverse - (2.16) 

See e.g. [11] for further details. This leads to 



j Sm » = * (2.17) 

where $| = is the South pole projection. 

The minus sign in (2.17) arises because stereographic coordinates are related by inversion, 
which has negative determinant. The expression for the strain tensor Ds m n is 

5mn (i + i? s 2 ) 2 ar m ax n 
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There is a sign ambiguity when defining the Jacobians Jn and Js- Here we have chosen 
the Jacobian such that the standard hedgehog has positive Jacobian. In fact, for B = 1 the 
hedgehog can be written as 

X ' 

a = cos(/(r)) and 7r» = — sin(/(r)), (2-18) 

r 

where r = ^x\ + x 2 , + x\ and / is a real radial shape function that can be determined 
numerically. Asymptotically, /(r) ~ n — as r — )■ with A > and /(r) ~ ^ as r — )■ oo 
with C ~ 2.16, see e.g. [12]. Near the origin, a ~ —1, so we project from the north pole 
and obtain det Jn ~ ^4 3 > 0, whereas as r — >■ oo, a ps 1 we project from the south pole and 
Js ~ > 0. These expression of the Jacobians (Jn and Jg) are useful when examining 
the behaviour near the north and south pole, which corresponds to the vacuum and the 
anti- vacuum respectively. 

For analysis of numerical data later in the paper, we can express the Jacobian directly in 
terms of the four pion fields as 

(J Nr /) = (2.19) 
ox m 

dn^ 7r^ da 
dx m (1 - a) dx m 
dirf 1 ^ tiV dix v 



S—^ VI ' 71 1/71 

dx m ^ a(l - a) dxj 



where we used the identity 

° dxj ^-/ 7r Qxn ' (2.22) 

which is derived by differentiating the normalisation condition. 
As a check, we can evaluate 

D mn = J m ^{x)J n (x)5^ l , (2.23) 

and obtain 

dn^ dn u / _ n fl n p 
dx m dx n V ^ a 2 

The term in the brackets corresponds to the induced metric in (a, TTj) coordinates. Using the 
identity (2.22) we obtain the more familiar expression 

da da d^ dir v m ok\ 

OX m OX n OX m OX n 



"--ttttV'-^- (2 - 24) 



With this expression, equation (2.9) gives the pion field version of the Skyrme Lagrangian 
which is often used for numerical simulations, see e.g. [32]. It can be shown by direct 
calculation that det (Jn) = — det (Jo) /a where 

<9tT,- 

J ° = &? < 2 ' 26 > 

Hence det J N vanishes when det J does for a ^ 0, therefore Jn and J have the same singular 
surfaces. 
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3 Rational map ansatz 



The rational map ansatz [9] makes use of the feature that a point in IR 3 can be written 
in polar coordinates (r, z) where the angular coordinate z is represented by a point on the 
Riemann sphere via z = e^tan |. The ansatz then takes the form 

U(r, z) = exp(if(r)n(z) ■ r), (3.27) 

where f(r) is a real valued profile function with the boundaries f(0) = it and f(oo) = 0; 
r = (n, T2, T3) is the triplet of Pauli matrices and n is the unit vector 



n(z) = 1 + TO | 2 W + R(*)AR{z) ~ AC*)), 1 " 1^)1 )• 

-R(z) is a holomorphic rational map between Riemann spheres, and is given by the two 
polynomials p(z), q(z), 

**) = 44- 

q{z) 

Then the baryon number is equal to the algebraic degree of R(z), 

deg (R(z)) = max {deg(p(z)), deg(g(z))} . 

Such an ansatz is a suspension and gives rise to an isomorphism between ^(S* 3 ) and ^(.S 2 ). It 
is the choice of R(z) which replicates the polyhedral shape of Skyrmions, and the symmetries 
of the accepted numerical solutions. The benefit of the rational map ansatz is that it gives 
rise to the following three simple eigenvalues of the strain tensor, 

Ai = -f(r) (3.28) 

^2 = A3 - 



sin / 1 + 


z 


2 


dR 


r 1 + 


R 


2 


dz 



This now gives the simple radial energy functional 



E = Att I ( fV + 2E(/ /2 + l)sin 2 f+l S -^- + 2m 2 r(l- cos f) ) dr, (3.20) 



where 



1 

47T 



1 + 


z 


2 


dR 


1 + 


R 


2 


dz 



Hdzdz 



1+ z 



2\2 ' 



(3.30) 



This energy functional can be easily minimised by choosing the correct degree polynomials 
p(z),q(z) which minimise X, then numerically minimise the profile function f(r). The min- 
imum energy solutions found using this method only exceed the non-symmetry numerical 
solutions by about 3%, see e.g. [HJ H31 EH]- 

This holomorphic ansatz is very successful at capturing the major features of the Skyrme 
model and is a very useful technique to give initial configurations which are close to the 
minimum energy solutions for numerical minimisation. This avoids a numerically expensive 
collisions which were previously used to create appropriate initial conditions [TS1 fT3"] . A 
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feature of the holomorphic ansatz is that it locally preserves orientation. Hence there cannot 
be regions of negative baryon density. This constrains the possible configurations. We can 
extend the ansatz to non-holomorphic R(z,z) as in [8]. This now allows points in IR 3 with 
negative baryon density. In this case the three eigenvalues of the strain tensor are slightly 
more complicated and are 

Ai = -/'(r), (3.31) 

A 2 = ^l(\R z \ + \R- z \) 
r 

A 3 = ^J-{\R t \-\R g \) 



1 + 


\z\ 


2 


1 + 1 


R 


2 



1 + 


\z\ 


2 


1 + 1 


R 


2 



where R z and R z are the derivatives of R with respect to z and z, repectively. It is now 
apparent that for non-holomorphic (or non-antiholomorphic) R the angular strains A 2 , A3 
are no longer isotropic, and only when R z (or R s ) equal zero do we regain the previous 
holomorphic ansatz. The energy functional for this more general ansatz is 

E = j- J (^f 2 r 2 + 2J(f 2 + 1) sin 2 / + X^^ + 2m 2 r(l - cos ffj dr, (3.32) 

where 

J- h J fa MM (3-33, 

/ 2\ 2 

~\2\2 ' 



X = -J l W - M ^-^ IrT ^ p , (3.34) 



and 

« = s/ W^/M fa-M (^) 2 ) (3.35) 

Again f(r) is a profile function which is a solution of the ODE 

(r 2 + 2 J sin 2 /)/" + 2rf + {Jf 2 -J- X^^) sin 2f - m 2 r sin / = 0, (3.36) 

with the boundary conditions /(0) = n and /(oo) = 0. 
As above, if we restrict to maps of the form 

R(z,z) P{zrz) 



q(z,z) 



where the z dependence is chosen in such a way as to preserve the symmetry of the original 
minimal energy rational map. Then the baryon number is generically equal to N1—N2, where 
iVi is the maximal holomorphic degree of (p, q) and N 2 is the maximal antiholomorphic degree 
of (p,q). 
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In (3.36) J replaces the role of B in the holomorphic rational map ansatz. This has a 
significant effect on the profile function /. To understand this we linearise (3.36) about r = 
and set f(r) = n — v{r) where u(r) <C 1. This gives the new linear ODE 



r 2 v" + 2rv — 2Jv + m 2 ru = 0, 



with the solution 



f(r) 



7T 



Cr^ 



-1+^1+8.7 ' 



(3.37) 



(3.38) 



For m = 0, equation (3.37) is a Cauchy-Euler equation, whereas for m ^ the solution is 



given in terms of a Bessel function. Numerical calculations show that as the pion mass is 
increased, J becomes larger and hence / decays slower about the origin. We can now play 
the same trick for the limit of r going to infinity. Here we set f(r) = e(r) where e(r) <C 1, 
giving the linearised equation 



2 // i o I 



r e 



2Je - m 2 re = 0. 



This gives the solutions 

f(r) 



1 + y/lW 

Cr V a— 



Ce" 2 ™^ 4(1+8J)-1 ^ 



r 3/4 



if m = 0, 
if m 0. 



The effect of the mass term is to make the solution more localised around the origin. 



4 The B = 3 Skyrmion 

It has been known for a long time that the baryon density for the B = 3 Skyrmion is tetrahe- 
drally symmetric [TOJ [15] . A lot of its features can be explained by a tetrahedrally symmetric 
rational map [9]. But, as shown in [8], the non-holomorphic rational map which allows 
negative baryon density has lower energy that the holomorphic map. The non-holomorphic 
rational map gives rise to four singular tubes (det(J) = 0) which start at the origin, pass 
through the faces of the tetrahedron then go off to infinity. It has been shown [8] that there 
are three folding lines equally spaced along these tubes. Also, in the centre of these tubes 
there are regions of negative baryon density. This inspired us to understand the form and 
amount of negative baryon density in actual minimum-energy Skyrme solutions. 

The family of rational maps for the B = 3 Skyrmion with the correct symmetry is 

R = p^se + P2 sine 

qi cos + q 2 sin 



where the polynomials Pi,p 2 , <?i and q 2 are 

Pi(z, z) = i\/3z 3 z + i\/3z 3 — zz — 1, Qi( z , z) = + z 3 — i\/3z 2 z — i\/3z, 



p 2 {z, z) = z" — 2i\/3z 2 + 1, Q2(z, z) = —z^z + 2i\/3z 2 z — 



v'l 



(4.40) 



With a simple numerical scanning algorithm we find that the family of rational maps in (4.39 ) 



attains its minimum energy for 8 = 0.154. This gives E/B = 1.161. Note this is lower than 
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Figure 1: Baryon number three minimum-energy solution. The red surface is a level set of 
constant positive baryon density. The yellow surface is a level set of constant negative baryon 
density. 



the minimum energy for the holomorphic ansatz (E/B = 1.184). This value of 9 is slightly 
lower than that found in [8]. This difference is believed to be due to numerical accuracy. 

To find the B = 3 minimum energy solution we numerically minimised this rational 
map ansatz on a lattice of 250 3 points, with spacing Ax = 0.08 using fourth order accurate 
derivatives. For m = we find E/B = 1.146 and the negative baryon density B_ = 4.5 x 10 -5 . 
This is in reasonable agreement with B- = 9.25 x 10~ 5 found for the non-holomorphic rational 
map for 9 = 0.154. A surface of constant baryon density (ie. constant Jacobian determinant) 
is displayed in figure [T] for this solution. 

To truly capture these regions of negative baryon density we were forced to have a very 
large box. It is known that asymptotically the B = 3 Skyrmion decays as a B = — 1 
Skyrmion [T7]. The B = 1 Skyrmion can be calculated from the highly symmetric rational 
map, R = z, known as the hedgehog ansatz. Then the equations of motion ( |3.36 ) can be 
linearised to obtain asymptotic behaviour of the radial profile function / ~ C/r 2 as r — > oo 
for some constant C. Substituting this into the radial energy-density (3.32) and baryon 
density (3.35), we find that the energy-density scales as 1/r 4 and the baryon-density as 1/r 7 . 
We have checked numerically that choosing U(x) = I2 on the boundary of the numerical box 
contains the baryon-density within the accuracy of our numerics. However, for an evaluation 
of the total energy to a high order of accuracy contributions from outside the box need to be 
included, see e.g. [18] for further details. 

Figure [T] shows the regions of negative baryon density of the minimum-energy baryon 
number three Skyrmion. A point of particular interest is the tetrahedron of negative baryon 
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> 



(a) Yellow is a level surface of det(J(a?)) 
the red is surface of the anti-vacuum a ~ 



(b) Preimages of the anti-vacuum. Red is where 

and sign( J(x)) = —1 and Ivory is where sign( J(x)) = 

1 1. 



Figure 2: Minimisation of the central region of the B = 3 non-holomorphic rational map. 



density in the centre of the Skyrmion. Each of the corners of this dual-tetrahedron gives rise 
to tubes of det(J) = 0, which pass though the faces of the tetrahedron with positive baryon 
density. However, the singular surface tubes may or may not be connected to the dual 
tetrahedron singular surface in the centre. With the aim of trying to understand how these 
singular surfaces are connected we performed a numerical simulation over a box centred at 
the origin with (Ax = 0.01, n = 180) sides about 9% of the large box. We chose the boundary 



conditions, on the surface of the box, to be the rational map (4.39) with the minimum-energy 
profile function found numerically which is a good approximation to the exact solution. Our 
intention is to present evidence to understand the singular structure about the origin of the 
B = 3 Skyrmion. The result of our numerical calculation is presented in figure [2j 

Figure 2a shows the folding tubes in yellow. Contrary to the conjecture in [8] the tubes 
do not seem to pinch off to singular points. This is evidence that there are no swallowtails 
in the B = 3 Skyrmion and will be discussed further in the next subsection. Figure 2b 
is also interesting. It shows that it is energetically favourable for the B = 3 Skyrmion to 
'create' two more pre-images of the anti- vacuum, one with positive orientation and one with 



negative orientation. The rational map ansatz (3.27) has the anti- vacuum as a suspension 
point. Hence for the rational map ansatz there are only three pre-images of the vacuum, all 
at the origin and with positive orientation. These extra pre-images are very similar to the 



monopole zeros of the Higgs field in [19] where they find five monopole zeros, four positive 
and one negative. As outlined in [Tj5], one would naively suspect that a positive orientation 
point would annihilate with a negative orientation point. However, the configuration may be 
stabilised by tetrahedral symmetry. 
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The 6 which yields the minimum-energy ratio- (b) E/B in (3.32) for the optimal value of 9 as a 



map ansatz (4.39). 
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function of m. 



3.049 




X in (3.34) as a function of m for the optimal (d) J in (3.34) as a function of m for the optimal 



Figure 3: The minimum-energy E/B and 9 as a function of m of (3.32) 



4.1 The effects of the pion mass on the folding structure 

It is well known that the Skyrme theory describes atomic nuclei better when a pion mass 
term is included, see pEJ EOJ, [21] for a discussion of m / and [5J for the implications of 
spinning Skyrmions. This inspired us to investigate how the surfaces of det(J(:r)) = vary 
as a function of m. We first examine this question using the non-holomorphic rational map 
ansatz and then check our results using full field minimisation. As a side effect, this analysis 
provides a good test of the effectiveness of the rational map ansatz. 



In order to find the minimum-energy rational map in (4.40) for given m we minimised 



the numerical integral of (3.32) with respect to 9 using a standard search algorithm. The 
resulting function 9(m) is displayed in figure 3a for values of m between and 1. The 
energy per baryon number is shown in figure 3b Note that E/B increases monotonically 
as m increases which agrees with the findings in [2D]. A point worth noting is that it is 
energetically favourable for X to decrease slightly, see figure 



3c 



and for J to increase as m 
increases, see figure [3d] This has an effect on the behaviour of the shape function around 
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the origin as can be seen from equation (3.38). 
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(a) The total negative baryon density as a function (b) Area of the negative baryon density tubes as 
of m. a function of m. 



Figure 4: Data arising from the minimum-energy B = 3 rational map (4.39). 



The radial integral of the baryon density in ( 3.35[ ) can be evaluated exactly using the 
boundary conditions of f(r). Using the minimum-energy 9{m) we then integrated the negative 
baryon density £>_ over S 2 to give the total negative baryon density B_ as a function of m 



as displayed in figure |4aj It is apparent that as the mass m increases it is more energetically 
efficient for the B = 3 Skyrmion to have more negative baryon density, B-. 

In the rational map ansatz, the negative baryon density arise as tubes emanating from 
the faces of the tetrahedral polyhedron, very similar to the yellow tubes in figure [TJ To 
understand how these tubes change as a function of mass, we numerically integrated the area 
of a negative baryon density tube over S 2 . This is shown in figure 4b Hence, the rational 



map ansatz predicts that the tubes of negative baryon density increase in size roughly linearly 
with m. 

To proceed we numerically minimised the rational maps for the values of the pion mass 
given above using a full field minimisation. The results are shown in figures [5j It is worth 
discussing the actual full field minimum-energy B = 3 solutions. One should first note that, 
as expected and conventionally understood, the hole in the level-set of positive det (J(x)) 
becomes smaller for increasing m. This is well understood as the potential term forces the 
Skyrmion field to reach the vacuum value exponentially. Also, as m increases the singular 
tubes become more pronounced. This is clearly seen in figure [5] where the singular tubes are 
much more defined for m = 1 than for m = 0.1. 



As shown in figure 6a, the energy per baryon for the full field minimisation is about 2% 
lower than that for the rational map ansatz, for all m. This is expected, because as previously 
shown the rational map ansatz forces all of the pre-images of the anti-vacuum to be at the 
origin. 



As can be seen in figure [2a| the Skyrme field has a fifth pre-image at the origin with negative 
orientation showing that there is also negative baryon density at the origin. This region of 
negative baryon density is inside a tetrahedron, which is dual to a surface of positive baryon 
density. In this central region a ~ — 1, this is where the potential energy (the pion mass 
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(e) m = 0.5 (f) m = 0.6 



term) is maximal hence as m increase the negative baryon density at the origin decreases. To 
verify this we numerically integrated over the central tetrahedron of negative baryon density, 



13 





(g) m = 0.7 




(h) m = Oi 




(i) m = 0.9 (j) m = 1 

Figure 5: B = 3 Skyrmion of varying m. Red is a level set of B = const > 0, yellow is B = 0. 



this is shown in figure [6cJ This shows that indeed as m increases the total negative baryon 
density in the central dual tetrahedron decreases. This gives a heuristic explanation why the 



total negative baryon density initially decreases, shown in figure [6b| as m starts to increase 
from m — 0, then after m ~ 0.7 the negative baryon density starts to increase again. This is 
reflected in the non-holomorphic rational map ansatz. This trend is further verified because 
we also found the minimum-energy B = 3, m = 2 solution. This solution was found to have 
E/B = 1.551,5_ = 4.1 x 1CT 5 and a £?_ dual tetrahedron total baryon density of 1.2 x 1CT 5 . 
Figure [6c] is not very smooth. This is most likely due to the numerical grid and the algorithm 
which identified the edge of the dual tetrahedron. 



4.2 Expansion around the origin 

The tetrahedral symmetry of the B = 3 Skyrme field poses stringent restrictions on the 
allowed terms in a Taylor expansion around the origin. Here we calculate the allowed poly- 
nomials and estimate the relevant coefficients from the numerical solution. 
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Figure 6: B = 3 full field minimum-energy solutions data. 



The rational map of the B = 3 Skyrmion is given by 

„. , y/3iz 2 — 1 , 

R(z) = ~ 3 7^7" (4-41) 

z 6 — y iiz 

which is Td tetrahedrally symmetric using the same orientation as in [5]. This symmetry is 
generated by a C 2 symmetry 

z i— > —z, R (-)■ -R, (4.42) 

and a C 3 symmetry 

iz + 1 n -iR + i ,. Anx 

flh+ TR+T' (443) 

for the tetrahedral symmetry T, together with an additional reflection symmetry 

z i—)- , i? H- = . (4.44) 
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The corresponding symmetries in component notation for the C 2 generator is 

(Xl,X 2 ,X 3 ) !->■ (-Xi, -X 2 ,X 3 ) , (7Ti,7r2,7r 3 ) f-^ (-7Ti, -7r 2 ,7r 3 ) . (4.45) 

and 

(xi, x 2 ,x 3 ) i-> (x 2 ,x 3 ,x 1 ) , (71-1,7^,^3) i-> (7r 3 ,7ri,7r 2 ) (4.46) 
for the C3 generator. The reflection symmetry is given by 

(xi, x 2 , x 3 ) ^ (a?!, -x 3 , -x 2 ) , (7Ti, 7r 2 , tt 3 ) h-> (tti, -tt 3 , -tt 2 ) (4.47) 

Hence the field tt can be expanded around the origin as 

yz \ I x 3 
7T = ai I y I +bi I X2; + c i 2/ 3 I + c 2 r 2 

XT/ / V Z 3 





yz 

+ di I xzy 2 I + rf 2 r 2 I xz ]+..., (4.48) 

xy 

where the Cartesian coordinates are now denoted as (x,y,z) and r 2 = x 2 + y 2 + z 2 . The 
tetrahedral symmetry can be augmented to spherical symmetry by setting b\ — 0, c\ = 
and c?i = c? 2 = 0, which corresponds to the Taylor expansion of the hedgehog ansatz around 
the origin. This provides a useful check that we have implemented the tetrahedral symmetry 
correctly. 

In order to compare the Taylor expansion of 7r to the numerical solution, we have to 



calculate the coefficients in (4.48). We used the following approach. Setting y — and z = 
in (4.48) gives a polynomial in x. We fitted this polynomial to our numerical data using a 
least square fit, and this allows us to calculate a x and C\ + c 2 . Similarly, setting y = x, and 
z = gives ax, b±, c% + 2c2 and 2cf 2 - Finally, setting y = x and z = x also gives an equation 
in d\. In order to improve the approximation close to the origin, we fitted to a polynomial 
of higher degree. By plotting numerical data and approximation we found that the values 
di = —0.41, b\ = —2.1, ci = —2.2, c 2 = 1.2, d\ = 0.76 and d 2 = 0.64 are a reasonable 
approximation for —0.5 < x,y,z < 0.5. Note that the errors in particular for the coefficients 
d\ and d 2 are rather large. We can now evaluate the Jacobian of the map. At the origin 
x = 0, the Jacobian is non-zero, namely, det(Jo) = a\ < 0. The value of det Jo(0) = —0.055 



calculated numerically from the exact solution matches the value for the expansion (4.48). 

Figure [7] shows a plot of the surface det J = inside a cube of length 0.5. This clearly 
looks very similar to the surface det J = arising from the numerical solution displayed in 
figure |2b| Using our expansion, we can check whether the singular surface pinches off at a 
point. By setting y = x, and z = x in the Jacobian Jo we can show that det Jo = for 
x = —0.09 but there is no positive solution within the box. We were careful to include terms 
up to fourth order because the normal form of the swallowtail includes a fourth order term. 
In summary, we have deduced the following folding structure for the B = 3 Skyrmion. There 
are four folding tubes through the faces of a tetrahedron. These tubes smoothly connect 
to the corners of a dual tetrahedron at the origin. The folding surface are not intersecting 
each other, so there are no swallowtail singuarities in the B = 3 Skyrmion, contrary to 
the conjecture in [8]. By symmetry, the cusp lines are expected to lie on the edges of the 
tetrahedron, and there is some numerical evidence. 
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Figure 7: Singular surface det J = calculated using the polynomial expansion around the 



origin in (4.48). 



5 Conclusion 



This paper was motivated by the results of [8], where the authors found regions of negative 
baryon density in the rational maps for the B = 3 minimum-energy Skyrmion. For small m 
these regions of negative baryon density are very small, but we have been able to numerically 
verify their existence. Also we have discovered a tetrahedron of negative baryon density at the 
origin, which is dual to a tetrahedron produced as a level-surface of positive constant baryon 
density. The singularities corresponding to surfaces of zero Baryon density form four tubes 
which smoothly join up at the dual tetrahedron at the origin. Contrary to the conjecture in 
[8], there are no swallowtail singularities in the B = 3 Skyrmion configuration. We have also 
found that for the B = 3 Skyrmion there are five pre-images of the anti-vacuum. Four with 
positive orientation, on the vertices of a tetrahedron, and one with negative orientation at 
the origin. This behaviour is also seen in monopoles [19|. 

The authors of [8] did not find any regions of negative baryon density in the rational map 
ansatz for the B = 4 minimum-energy solution. This has been verified here numerically. 
Furthermore, assuming octahedral symmetry we have shown in appendix that there are no 
regions of negative baryon density around the origin. These results are also consistent with 
the instanton ansatz [22]. 

It has already been discussed that regions of negative Jacobian-determinant occur for 
the charge three instanton |22j. This is signifiant because there is a BPS extended Skyrme 
model which can be derived from Yang- Mills instantons [23], which must contain regions 
of negative baryon density. This extended model has an infinite number of vector mesons 
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and rho mesons. There has also been research into a truncated version of this model [24] , 
where only a few extra terms are included. Understanding the form and distribution of the 
negative baryon density in these models would be very interesting. It should be noted that 
there is another BPS Skyrme model [25J, where the Bogomolny equation shows that the 
baryon density is proportional to the square root of the potential. Hence, if the potential is 
positive definite through space, so is the baryon density R 
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A B = A 



It has been shown [8] that even when the B = 4 rational map ansatz is extended to be 
non-holomorphic no negative baryon density is found. This is also seen in our full field 
minimisation of the B = 4 Skyrmion for m = and m = 1, see figure 8a and 8b, respectively. 
For m = we found E/B = 1.12, and for m = 1 we found E/B = 1.30. In both cases, up to 
numerical accuracy, we did not find any negative baryon density. 




(a) m = 



(b) m = 1 



Figure 8: minimum-energy full field solutions of the B = 4. Red is a surface of constant 
positive det(J). 



In the following, we show that the B = 4 minimum-energy octahedrally symmetric 
Skyrmion does not have regions of negative baryon density near the origin. This is where 
the non-holomorphic rational map ansatz is not a good approximation of the exact solution. 

1 This was pointed out by M. Speight. 
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As a starting point, we consider the rational map of the B = 4 Skyrmion is given by 

— 2-\/3zz + 1 

which is octahedrally symmetric. This symmetry is generated by a C4 symmetry: 

z h-> iz, -R — , (1.50) 
it 

and a C3 symmetry 

2^ — , R^e~ — R. 1.51 

-zz + l' V ; 

The corresponding symmetries in component notation for the C4 generator is 

(x 1 ,x 2 ,x 3 ) f-> (x 2 , -x 1 ,x 3 ) , (7ri,7ra,7r 3 ) •-»> (7ri,-7r 2 ,-7r 3 ). (1.52) 

and 

(x 1 ,x 2 ,x 3 ) H> (x 2 ,x 3 ,xi) , (7ri,7r 2 ,7r 3 ) i-> ^-^7Ti - ^\/37r 2 , ^\/37ri - ^7r 2 , 7T 3 J (1.53) 

for the C 3 generator. Hence the field "ir can be expanded around the origin as 

x 2 + y 2 - 2z 2 
n — bi \ -V3x 2 + V3y 2 


x 4 + y A — 2z 4 \ / 2x 2 y 2 — x 2 ,2 2 — y 2 ^ 2 

+ di | -v^a; 4 + v^Z/ 4 j + d 2 I ~V3x 2 z 2 + y/3y 2 z 2 | + 

A similar expansion as r tends to infinity has been performed in [TS]. We can now evaluate 
the Jacobian of the map. At the origin x = 0, the Jacobian is identically zero. Hence the 
singularity is clearly non-generic since for a generic singularity the Jacobian has rank 2. The 
determinant of the Jacobian can be expanded into O— symmetric polynomials. The lowest 
order term is 

det(J) = V3b 2 Cl (x 2 y 2 + y 2 z 2 + 2V) + . . . 

Hence the singularities are on the three coordinate axis which meet at the origin. Further- 
more, sufficiently close to the origin det(J) is positive, and hence there is no negative baryon 
density. 
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